Gamow-Teller strength for the analog transitions to the first T=l/2, J'^ = 3/2 states 
in ^"^C and ^'^N and the implications for Type la supernovae. 
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The Gamow-Teller strength for the transition from the ground state of ^^C to the T=l/2, J'^ = 
3/2" excited state at 3.51 MeV in is extracted via the ^^C(^He,t) reaction at 420 MeV. In 
contrast to results from earlier (p,n) studies on ^"^C, a good agreement with shell- model calculations 
and the empirical unit cross section systematics from other nuclei is found. The results are used to 
study the analog ^^N(e~, !/e)^^C reaction, which plays a role in the pre-explosion convective phase 
of type la supernovae. Although the differences between the results from the {^Ke,t) and {p,n) data 
significantly affect the deduced electron-capture rate and the net heat-deposition in the star due to 
this transition, the overall effect on the pre-explosive evolution is small. 

PACS numbers: 21.60.Cs, 24.50.+g, 25.40.Kv, 25.55.e, 25.60.Lg, 26.50+x, 27.20.+n 



I. INTRODUCTION 



Type la supernovae (SNe la) are thought to be ther- 
monuclear incinerations of carbon-oxygen white dwarf 
stars after accretion from a companion has increased their 
mass to near the Chandrasekar limit (for a review, see e.g. 
[1]). It is unclear, however, what type of binary systems 
produce SNe la and how the structure of the progenitor 
affects the outcome of the explosion. Observational stud- 
ies show that the peak brightness of SNe la, and hence 
the mass of ^^Ni synthesized in the explosion, are cor- 
related with the properties of the host galaxy [2]. One 
proposed explanation for this is that the ratio of ^^Ni 
to stable nickel and iron isotopes depends on the elec- 
tron abundance Ye in the white dwarf at the time of the 
explosion. A reduction in shifts the isotopic distri- 
bution of iron-peak ejecta to more neutron-rich isotopes 
[3] and decreases the mass of '^^Ni ejected, which in turn 
reduces the peak brightness of the supernova. In the pre- 
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explosion white dwarf, is set by the abundance of trace 
nuclides such as ^^Ne [4]. 

In recent works [5, 6] it was pointed out that reactions 
during the pre-explosion "simmering" phase (during the 
~ lO^yr prior to the explosion) also reduce Ye- During 
this period, protons formed from ^^C{^^C,p)^^Na, pre- 
dominantly capture onto ^^C to form /3-unstable ^^N. 
The ^3N(e~, i^e)^^C reaction thus plays a role. Since the 
stellar density is p = (1~3) x 10^, the electron Fermi en- 
ergy is sufficiently high (« 5 MeV) to allow electron cap- 
ture on the ground state (spin-parity — l/2~) 
via an allowed Gamow-Teller (GT) transition into the 
J'^ = 3/2- state of at 3.68 MeV [5], in addition to 
the capture to the ^^C ground state (J'^ = 1/2"). Note 
that captures into states at higher excitation energies are 
not important for the range of densities found in white 
dwarf stars, so the heating from this reaction is set by 
captures into this excited state and the ground state. As 
discussed briefly in Ref. [5] , where a shell model calcula- 
tion was used to estimate the contribution of the excited 
state, the increased capture rate does not strongly effect 
the reduction of Y^, which is instead governed by the 
production of ^^N, not by the faster (under these con- 
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editions) electron capture. The capture into the excited 
state does contribute to the effective heat per ^^C + ^^C 
reaction during simmering, however, and so a more accu- 
rate determination of the branching into the excited state 
is useful for future studies of the simmering stage of SNe 
la, including those focusing on the hydrodynamics. Hy- 
drodynamical simulations of the simmering phase [7, 8] 
require an accurate determination of the effective heat 
released per ^^C +^'^ C reaction, since this heating sets 
the amount of ^^C consumed, and hence the neutroniza- 
tion of the white dwarf, prior to the explosion. Motivated 
by the previous work [5], this paper incorporates exper- 
imental data into the determination of the electron cap- 
ture rate into the excited state and the related increase 
in heat deposition. 

The weak transition strength from the ^■^N ground 
state to the ^'^C ground state is well determined using 
the experimental log/t value of 3.6648 ± 0.0005 for (3'^- 
decay [9]. Because the excitation energy of the 3.68 
MeV state is higher than the Q-value for /S^-decay (2.22 
MeV) , no such data are available for the transition to this 
state. Therefore, in Ref. [5], the GT transition strength 
from the ^'^N ground state to the excited state at 3.68 
MeV in was calculated (B(GT)=1.5) using the shell- 
model code OXBASH [10] employing the CKII interac- 
tion [11] in the p-shell model space, taking into account 
a phenomenological quenching factor of 0.67 [12] for the 
Gamow- Teller strength ^. 

Experimental information about the Gamow- Teller 
strength can be extracted, however, by studying the ana- 
log transition from the ^^C ground state to the — 
3/2~ state at 3.51 MeV in -'^^N, using a (p,n)-type charge- 
exchange reaction and assuming isospin symmetry (the 
validity of this assumption is discussed below). In the 
simplified level-scheme of relevant A=13 nuclei depicted 
in Fig. 1, these analog transitions are labeled 5 and 6. 

The charge-exchange reaction studies rely on the pro- 
portionality between B(GT) and the charge- exchange 
cross section at zero momentum transfer (g = 0) (see 
details below). Several experiments on ^'^C have been 
performed using the (j),n) reaction at beam energies be- 
tween 120 MeV and 200 MeV. In the earlier experiments 
[13-15], aB(GT) of 0.82 ±0.05 ^ was extracted. The un- 
usually large discrepancy between the experimental value 
and the theoretical one (B(GT)=1.5, taking into account 
quenching) was the topic of considerable debate [14-16]. 
More recently, a measurement of the ^^C{p,n) reaction 
at 197 MeV gave B(GT)=1.06±0.05 [17]. This value is 
closer but still significantly lower than the shell-model 
predictions. Because of the significance of this transition 
for the calculations performed for SNe la, further inves- 
tigation is warranted. 



^ This value is calculated by squaring the quenching factor of the 
free-nucleon operator given by 5s = 1 — 0.19( jg)'^ '^^ [12]. 

2 This is the value from Ref. [15]. The B(GT) values in Refs. 
[13, 14] are within the uncertainties consistent with this value 
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FIG. 1: Isospin analogous transitions in the A=13 and 
Tz = ±1/2, ±3/2 isobar system are schematically shown. The 
Coulomb displacement energies are set to zero so that the 
isospin symmetry of the system becomes clearer. For each 
state, spin-parity and isospin T are indicated. Excitation 
energies are given, where relevant for the present work, rela- 
tive to the ground state of the nuclei. Transitions labeled by 
1-4 are analogs and transitions with labels 5 and 6 are analogs. 
For transitions labeled with the log/t's are known from 
/3-decay experiments. 



As an alternative to the ^^C{p,n) reaction, we studied 
the ^^C(^He,t) reaction at 420 MeV and combined the re- 
sults with an earlier ^'^C('^He,t) experiment performed at 
450 MeV [18]. The data at 420 MeV were taken as part of 
a program to study the extraction of GT strength using 
the ('^He,i) probe [19] over a wide mass range. Since the 
targets used contained natural carbon (1.07% ^^C) only 
transitions to low-lying states in ^'^N could be studied. 
The transition to the T = 3/2, ,r = 3/2" state at 15.1 
MeV in ^"^N is important for strength calibration pur- 
poses. However, in targets made of natural carbon this 
state is inseparable from the strong ^^C('^He,i)^^N(g.s.) 
reaction since the Q-values are nearly identical. There- 
fore, the data from the older experiment, in which a 
highly-enriched ^'^C target was used, was needed as well. 
Comparing the two sets also provides a good way to eval- 
uate systematic errors in the extraction of cross sections. 

This paper is structured as follows. After briefly de- 
scribing the method of extracting GT strengths from 
charge-exchange data and the particular methods used 
for the A=13 analog transitions, the findings from the 
^^C('^He,i) experiments are presented and compared with 
the shell-model calculations. It is shown that the results 
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from the ^■^C(''Ho,t) experiments and the shell-model cal- 
culations arc consistent. We then briefly discuss the com- 
parison with the (p,n) results. Finally, we estimate the 
uncertainties in the electron-capture rates on ^'^N in the 
stellar environment and discuss the ramifications for the 
simmering stage of SNe la. 



II. EXTRACTION OF WEAK TRANSITION 
STRENGTH FROM THE "C(^He,t) REACTION. 

For both the (p,n) and (''Hc.f) reactions, the extraction 
of B(GT) values from the data is based on their propor- 
tionality to the cross section at zero momentum trans- 
fer {q = 0) at sufficiently high beam energies (> 100 
MeV/nucleon). In eikonal approximation, this propor- 
tionality is written as [15]: 

^(g = 0) = KN^r\Jar\^BiGT) = a,rB{GT). (1) 



Here, K 



where -E'i(/) is the reduced energy in 



the incoming (outgoing) channel, N^t is the distortion 
factor defined by the ratio of the distorted-wave to the 
plane-wave cross sections: | J^tI is the volume-integral of 
the central (tt interaction. The factor KN\J^t\'^ is re- 
ferred to as the unit cross section, a^T- For excitations 
from a particular target nucleus, the unit cross section is 
usually determined experimentally, using transitions for 
which the B(GT) is empirically known from /3-dccay. 

The proportionality of Eq. 1 is not perfect. In partic- 
ular, contributions from incoherent and coherent AL=2, 
AS=1 to GT transitions lead to deviations. The inco- 
herent contribution, associated with a total angular mo- 
mentum transfer to the relative motion between the tar- 
get and projectile AJ/j= 2, can be extracted from the 
data, since its angular distribution peaks at finite angle 
and has a minimum at 0°. Usually, the AJij = 2 con- 
tribution to the 0° cross section is no more than a few 
percent of the total cross section. The coherent contri- 
bution, largely due to the non-central tensor-interaction 
and the largest source that breaks the proportionality of 
Eq. 1 [20-22] , cannot easily be determined from the data 
since it does not strongly affect the angular distribution 
at forward scattering angles. Its effect on the cross sec- 
tion must be estimated by theory. Such estimates have 
shown to provide reasonable estimates for the breaking 
of proportionality for (■^IIe,i) reactions on ^^Mg [20] and 
58Ni [21]. 

For Fermi transitions, a relationship similar to Eq. 1 
exists between j^{q = 0) and the Fermi strength B(F) 
with unit cross section (7^, but with a distortion factor 
Nt and volume integral of the central r interaction ] J^j- 
Since no spin transfer is involved, breaking of the pro- 
portionality due to AL=2 contributions is not an issue. 

Here, we are particularly interested in the excitation 
of the 3/2~ state at 3.51 MeV (transition 6 in Fig. 
1). The imit cross section used to convert the cross 
section at g = to strength can be calibrated in two 



ways. The first method relies on the use of the tran- 
sition to the -'^^N ground state, for which the strength 
is known from /3+ decay data [9]. This transition con- 
tains both a Fermi component (non-spin-flip, AS=0) and 
GT component (AS=1). Under the assumption that the 
Fermi strength B{F) exhausts the full Fermi sum rule 
{S-{F) - S+{F) = {N - Z) = 1)"^, the strength of the 
GT component B{GT) = 0.207 ± 0.002 can be found 
using: 



B{F) 
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with (|A) = 1.264 ± 0.002 [24] and K/g^ = 6147 ± 7 
[25]. However, to extract the unit cross section for the 
GT component only, the unit cross section for the Fermi 
component must be known with good accuracy, which 
makes this transition harder to use as a calibration tool. 

Fortunately, a second method to independently cali- 
brate the GT unit cross section is available, namely us- 
ing the transition to the '^^lSl{3/2- ,T = 3/2) state at 
15.1 MeV (labeled 3 in Fig. 1). This transition has 
three analogs, labeled 1,2 and 4 in Fig. 1. For transi- 
tion 1 (i3B(/3-)i3c) and transition 4 (1-''0(/3-)13n) the 
lifetimes, and thus B(GT) values, arc known from /?- 
decay experiments. After correcting for the difference 
in spin and isospin factors ^, the B(GT) for transition 3 
is estimated to be 0.213 ± 0.001 using the B(GT) from 
transition 4 and 0.237 ± 0.001 using the B(GT) from 
transition 1. The mirror asymmetry in the GT /3-decay 
matrix elements can be attributed to isospin breaking 
due to the Coulomb interaction. In Ref. [26], the ra- 
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M^C-'OC^/2-,g.s.)^^''N{l/2~,g.s.)) 



was estimated 



M2(i3B(;V2-,3.s.)^l'"C'(l/2-,3.s.)) 

to be 0.925 for the radial overlap (Method I of Table 
1 from Ref. [26]). Combined with value for isospin 
mixing in the p-shell for this transition of 0.992, R be- 
comes 0.918, which is close to the measured asymmetry 
of [;;2i3±o.ooi ^ o.899±0.006. In Ref. [15], the B(GT) for 
the transition 3 of Fig. 1 was set to 0.23±0.01, where the 
error was chosen to represent the minor ambiguity due 
to isospin symmetry breaking and, for the sake of consis- 
tency, we used the same number in the current analysis. 

We used the methods of Ref. [26] to estimate the 
isospin symmetry breaking between transitions 5 and 6 
in Fig. 1. The calculations are based on the Skyrme 
Hartree-Fock radial wave functions obtained with the 
Skx interaction [27]. The p-shell isospin breaking was 
obtained in proton-neutron formalism with the charge- 
dependent interaction of Ormand and Brown [28]. The 

rr>«i,lt fnr 7? - ("Af(l/2- ,g.^.)->i^C(3/2- ,3.68 MeV)) _ 
ICbUlL iOl n — M2(i3c(l/2-,g.s.)^i3Af(3/2-,3.51 MoV)) ~ 

1.043 for the radial overlap and 1.014 for the p-shell 



^ This assumption is valid on the level of 0.15(9)% [23] 

^ B(GT) is calculated from tho matrix element M(GT) via 

2Tf+l 



B(GT) : 



2(2Ji + l) " 



-M'^{GT) 
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isospin mixing for this transition for a total ratio R = 
1.06. Since this estimated asymmetry is comparable in 
magnitude to the uncertainties in extracting the B(GT) 
value for transition 6 via the charge-exchange reactions, 
isospin symmetry is assumed in the remainder of this pa- 
per. 

A minor disadvantage of using transition 3 to calibrate 
the GT unit cross section is that the Q-value (17.248 
MeV), and thus the minimum momentum transfer (at 
0°) for this transition are relatively large (0.18 fm~^). 
Since the cross section at g = must be used in Eq. 1 an 
extrapolation is required which is estimated theoretically 
(see below). However, the error in this extrapolation is 
small compared to the error in determining the GT unit 
cross section from the ground-state transition where the 
Fermi component has to be subtracted. 

After the determination of the GT unit cross section 
using transition 3, the B(GT) for transition 6 can be 
extracted. Moreover, since the B(GT) for the ground- 
state transition is known, the cross section for the GT 
component can be estimated and subtracted from the 
total cross section for the ground-state excitation. The 
remaining Fermi cross section is then used to determine 
the Fermi unit cross section. 



III. EXPERIMENTAL RESULTS. 

The i3c(''Hc,t) data were obtained at RCNP using the 
spectrometer Grand Raiden [29]. A 3-pnA, 420 MeV ^He 
beam bombarded a 4.1 mg/cm^ natQ target and a 1.0 
mg/cm^ ""*C8H804 (Mylar) target. Tritons were mea- 
sured in the focal plane of the spectrometer and mo- 
menta and angles reconstructed through a ray-tracing 
procedure. The spectrometer was tuned to run in the 
"off-focus" mode [30] of operation, so that the angles in 
both dispersive and non-dispersive direction could be ex- 
tracted with good accuracy. Center-of-mass scattering 
angles between 0° and 3.0° were covered with a resolu- 
tion of 0.2° (FWHM). The energy resolution was - 110 
keV with the thin Mylar target and ~ 150 keV with the 
""*C target due to larger differential energy loss of the 
"^He and triton in the thicker target. "^He beam particles 
were collected for the purpose of cross section calcula- 
tion in a Faraday cup placed in the inner bend of the 
first dipole magnet of Grand Raiden. Since the cross sec- 
tions for the two different targets were consistent within 
statistical errors and the total number of events collected 
a factor of ten larger with the ""*C target, the results re- 
ported here are those taken with the ""*C target. In the 
same beam-time period as thcisc; data were taken, data 
for other targets were obtained to establish an empirical 
relationship for unit cross sections as a function of mass 
number. The results of that study were presented in Ref. 
[19]. 

In Fig. 2a the excitation-energy spectrum of ^^N taken 

with the ""^'C is shown for center-of-mass scattering an- 
gles between 0° and 3°. Besides the ground state and 
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FIG. 2: a) Excitation energy spectrum of N taken with the 
Ciie,t) reaction at 420 MeV on a ""'C target. The ground 
state and excited state at 3.51 MeV are indicated, b) Dif- 
ferential cross sections for the ^^C(^He,t)"N(3.51 MeV) re- 
action from the present data (error bars are statistical only) 
and from Ref. [18]. The latter have been scaled by a factor 
0.88. Also included is the result from the DWBA calculation, 
for which the scale was adjusted to fit the data, c) Idem, 
but for the ^^C(^IIe,t)^^N(g.s.) transition. The contribution 
from GT and Fermi components to this transition is described 
in the text, d) Idem, but for the "C(^He,f)"N(15.1 MeV) 
transition. Only data from Ref. [18] is available. 
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the state at 3.51 MeV, the 1/2+ state at 2.36 MeV can 
also be seen. Fig. 2b shows the differential cross sections 
for the transition to the 3.51 MeV state. The data from 
Ref. [18] are also shown after scaling by a factor of 0.88. 
The effect of the slight difference in beam energy (420 
MeV for the present data and 450 MeV for Ref. [18]) on 
the absolute cross section was estimated in DWBA (de- 
tailed below) and was found responsible for about half of 
the required adjustment. The need for further scaling is 
most likely due to systematic errors in the beam integra- 
tions and target thicknesses. ^. Note that another state 
(J'^ = 5/2+) exists at an excitation energy of 3.55 MeV, 
which cannot be separated from the state at 3.51 MeV. 
The transition to the 3.55 MeV state is of dipole nature. 
If it had a significant cross section it would have dras- 
tically changed the angular distribution shown in Fig. 
2b, since dipole transitions have a minimum and maxi- 
mum at scattering angles of 0° and ~ 1.2°, respectively. 
We, therefore, concluded that the contribution from this 
transition can be ignored. 

Fig. 2c shows the measured differential cross sections 
to the ^^N ground state from both the present data and 
from Ref. [18]. The latter has been scaled by the same 
factor 0.88 described above. A slight discrepancy can 
be seen between the data sets, which is likely caused by 
a reduction of the relative contribution from the Fermi 
component to the excitation of this state at the higher 
beam energy. The strength of the central r interaction 
drops significantly as a function of beam energy [31, 32] 
and is apparently noticeable even when the beam energy 
changes only by 7%, as is the case for the comparison 
between the data sets at 420 MeV and 450 MeV. 

Finally, in Fig. 2d, the differential cross section for the 
transition to the -'^'^N state at 15.1 MeV is shown after 
scaling by the factor 0.88. Only cross sections from Ref. 
[18] are available in this case because of the contamina- 
tion of ^"^G in the targets used for the later experiment 
(see above). It should be noted that the ^^C target used 
in Ref. [18] was enriched to at least 99% and that the 
contamination from the -"^^N groimd state into the ^''N 
excited state at 15.1 MeV is no more than 4% based on 
the cross section measured for the transition to the ^-^N 
ground state in the data obtained with the ""*C target. 

To determine the cross section at g = 0, first the cross 
section at 0° was obtained from the measured differen- 
tial cross sections shown in Figs. 2b-d. For the excita- 
tion of the states at 3.51 MeV and 15.1 MeV, this was 
done by fitting the differential cross sections calculated in 
distorted-wave Born approximation (DWBA) to the data 
using a single scaling factor. The DWBA cross section at 
0° multiplied by the fitted scaling factor was then used 



Although the accuracies of beam integration and target thick- 
ness are not necessarily better in the newer data than in Ref. 
[18], we chose to scale the older cross sections to the newer ones 
to maintain consistency with Ref. [19], rather than talcing the 
average of the two data sets. 



as the estimate for the experimental cross section at 0° 
and its uncertainty determined from the error in the fit. 

The DWBA calculations were performed using the 
code FOLD [33]. in this code, the Love-Franey nucleon- 
nucleon interaction [31, 32] is double-folded over the 
projectile-ejectile and target-residue transition densities. 
Because the ^He and triton are composite particles, a 
short-range approximation as described in Ref. [31] was 
used for the exchange terms in the potential. For ^He 
and '^H, densities used in the folding were obtained from 
Variational Monte-Carlo results [34]. One-body transi- 
tion densities (OBTDs) for the transitions from -'^^C to 
were calculated with the code OXBASH [10] employ- 
ing the CKII interaction [11] in the p shell- model space. 
Radial wave functions were calculated using a Woods- 
Saxon potential. Binding energies of the particles were 
determined in OXBASH [10] using the Skyrme SK20 in- 
teraction [27]. It should be noted that the B(GT) values 
for the ground, 3.51-MeV and 15.1-MeV states, calcu- 
lated with the CKII interaction, were 0.190, 1.50 and 
0.228, respectively, after correcting for the quenching fac- 
tor of 0.67. These are slightly different, but a little closer 
to the experimental values from /3-decay for the ground- 
state and 15.1 MeV excitations, from those calculated 
with the WBT interaction [35] in Ref. [17] (B(GT) val- 
ues of 0.1746, 1.3444 and 0.2849, respectively). 

Optical potential parameters extracted from ^He elas- 
tic scattering on ^^C [18] were used in the DWBA cal- 
culation. Following Ref. [36], the depths of the triton 
potentials were calculated by multiplying the depths of 
the ^He potentials by 0.85, while leaving radii and dif- 
fusenesses constant. 

For GT transitions, in principle a fit in which the AL = 
and incoherent AL = 2 pieces are independently scaled 
is required. However, inclusion of an incoherent AL = 
2 component changed the contribution from the AL = 
component to the forward-angle cross section by less 
than 1%. Since this change is smaller than uncertainties 
due to statistical error margins, the incoherent AL = 2 
contributions are not shown in Figs. 2b-d. Such small 
contributions from the incoherent AL = 2 components 
are consistent with the DWBA calculations. 

For the transition to the ground state the procedure 
is complicated by the contribution from both Fermi and 
GT components whose angular distributions at forward 
angles are very similar. The contribution from the GT 
component must, therefore, first be fixed using the unit 
cross section determined from the transition to the ex- 
cited state at 15.1 MeV. 

In Table I the results from the extraction procedure 
are summarized. The B(GT) values derived from the 
log ft values obtained from /3-decay are given in row a), 
where available. For comparison, the B(GT) values cal- 
culated in the shell- model arc shown in row b). These 
strengths have been multiplied with the GT quenching 
factor of 0.67 [12]. A good consistency with the values 
extracted from /3-dccay data is found for the transitions 
to the ground state and the excited state at 15.1 MeV. 
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In row c) the extracted differential cross sections at 0° 
obtained from the data as described above are displayed. 
The cross sections at g = are shown in row d). These 
were obtained from the cross sections at 0° by multiply- 
ing them with the ratio ^{q = 0)/;^(0°) . The differ- 

ential cross sections in the numerator and the denomina- 
tor are both calculated in DWBA. These ratios are 1.01, 
1.03 and 1.25 for the transitions to the ground state (the 
ratio is the same for the Fermi and GT components), the 
excited state at 3.51 MeV and the excited state at 15.1 
MeV, respectively. The GT unit cross section (shown in 
row e) is calculated by dividing the cross section at g = 
in row d) for the excited state at 15.1 MeV by the cor- 
responding B(GT) shown in row a). This GT unit cross 
section is then used to calculate the B(GT) for the ex- 
cited state at 3.51 McV, given in row f). This B(GT) is 
close to the corresponding shell-model value provided in 
row b). The error bar in the experimental B(GT) value 
is due to statistical uncertainties in the data, uncertain- 
ties in the fitting procedure to the DWBA calculations 
and also includes the uncertainty in the unit cross section 
derived from the 15.1 MeV state as quoted in row e). 

Finally, the Fermi unit cross section was derived from 
the transition to the ground state. To do so, the esti- 
mated GT contribution in row g) to this excitation must 
be subtracted first and this is done by multiplying the 
GT unit cross in row e) with the B(GT) extracted from 
/3-decay (given in row a) for the ground-state transition. 
The remaining cross section, placed in row h) was then 
attributed to the Fermi transition. Since B(F)=1, this 
cross section is also the Fermi unit cross section, as shown 
in row i). 

As mentioned above, the proportionality between 

B(GT) and differential cross section at g = can be bro- 
ken if the interference between the AL = amplitude and 
AL = 2 amplitude due to the non-central tensor-r inter- 
action is strong. This has been noted before in studies 
of GT strengths of relevance for astrophysical purposes 
(see e.g. Refs. [21, 37]). Here, we follow the procedure 
previously used in Refs. [20, 21] to quantify the effects 
of such interference. Namely, by comparing the differ- 
ential cross sections calculated in DWBA with the full 
Love-Franey interaction with those in which the tensor- 
r component of the Love-Franey interaction is removed. 
The relative changes between these two values were com- 
parable to statistical error bars for the transitions to the 
excited states at 3.51 and 15.1 MeV, but large for the GT 
component of the transition to the ground state; includ- 
ing the tensor-T interaction increases the cross section by 
15%. A corrected estimate for the GT contribution to the 
cross section of the transition to the ground state, which 
takes this increase into account, is given between brack- 
ets in row g) of Table I. Consequently, the estimate for 
the component in the cross section for the transition to 
the ground state due to the Fermi component decreases 
to the values shown between brackets in rows h). As a 
result, the Fermi unit cross section also decreases to the 



value shown between brackets in row i). 

The extracted GT unit cross section of 20 ± 1 mb/sr 
agrees well with the empirical relationship (cto-t = 
109A^'^ '^''' = 20.6 mb/sr for A=13) of unit cross section 
as a function of mass number [19]. The empirical rela- 
tionship for the Fermi unit cross section (ar = 72A~^ °^) 
gives 4.75 mb/sr for A=13. This value is consistent with 
the extracted Fermi cross section if the correction to the 
GT component for the transition to the ground state due 
to the interference with the tensor-r interaction is applied 
(5.1 ± 0.5 mb/sr). Together with the previous studies in 
Refs. [20, 21], this provides further indication that the 
procedure for estimating the proportionality breaking de- 
scribed above is reasonable. 

The unit cross sections can also be calculated from the- 
ory, by dividing the GT and Fermi cross sections calcu- 
lated in DWBA by the shell-model B(GT) and B(F)=1, 
respectively. Values of (7^^ (theory) = 34.7 mb/sr and 
(Tt- (theory) = 6.31 mb/sr are found, which are both sig- 
nificantly higher than the experimental values. Such an 
overestimate of the DWBA cross section is seen over a 
wide target-mass range [19]. Although the reasons for 
this discrepancy are not well understood, a likely cause 
is the approximate treatment of the exchange contribu- 
tions which has been shown to lead to an overestimation 
of the cross sections in the DWBA calculations [38, 39]. 
Other systematic discrepancies could be due to ambigu- 
ities in the optical potential parameters employed in the 
DWBA calculation and possible density-dependences of 
the effective nucleon-nucleon interaction. 



IV. COMPARISON WITH "C(p,n) DATA 

In ^•^C(p,n) experiments detailed in Refs. [13-15], the 
GT unit cross section was determined using measure- 
ments of the ground-state cross section, combined with 
a measurement of the spin-flip probability Snn- When 
the extracted unit cross section was used to determine 
the B(GT) for the transition to the excited state at 
3.51 MeV, the afore-mentioned large discrepancy with 
the shell-model calculations was found. It was further- 
more noted that the GT unit cross section derived from 
the transition to the 15.1 MeV state was consistent with 
that determined using the ground-state transition and 
that the Fermi unit cross section was exceptionally large. 

In contrast to this, in Ref. [16], it was suggested 
that the proportionality between the cross section and 
GT strength was broken for the GT transition to the 
^^N(g.s.). An effective (p,n) operator was introduced in 
the shell-model calculations based on the description of 
quenching of GT strength in the sd-shell [40] . This mod- 
ified operator, designed to take into account medium ef- 
fects, contained a term (5p(87r)i/2[y(2) ^ g]A,/=i^^^ rj^^j^ 

term mediates contributions from AL = 2 amplitudes 
that interfere with the Ai = GT amplitudes. It was 

speculated [16] that the apparent enhancement of 5p in 



transitions going from j = I 



\ to 3 



\ (such as 
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TABLE I: Summary of the analysis for the ground state and excited states at 3.51 MeV and 15.1 MeV in N excited via 

the ^^C(^He,f) reaction. Note that the error margins in the experimental differential cross sections are due to statistical 
uncertainties in the fitting of measured angular distributions only. Detailed descriptions of each of the rows, labeled a)-i) in 
the table, are given in the text. 









state in ^^N 








1/2- g.s. 


3/2- 3.51 MeV 


3/2- 15.1 MeV 


a) 
b) 
c) 
d) 
e) 
f) 
g) 
h) 
i) 


B(GT)^_decay 

B(GT)cKii X 0.67" 

[ln(9 = 0l., (mb/sr) 

(Tctt (mb/sr) 
B(GT) 

[|n(9 = 0)]^s=i (mb/sr) 

[^(9 = 0)1^5=0 ("ib/s r) 
(Tr (mb/sr) 


0.207 ± 0.002 
0.19 

10.0 ± 0.5 

10.1 ±0.5 

4.2±0.2''(5.0±0.2') 
5.9±0.5'(5.1 ±0.5") 
5.9±0.5''(5.1 ±0.5") 


1.50 
26.5 ± 0.3 
27.3 ±0.3 

1.37 ±0.07 


0.23 ± 0.01 
0.23 
3.6 ±0.1 
4.5 ±0.1 
20 ± 1 



"The factor 0.67 represents the phenomenological quenching fac- 
tor for GT strength. 

''Without taking into account corrections due to the tcnsor-r in- 
teraction (for details, see text). 

'^After taking into account corrections due to the tensor-r inter- 
action (for details, see text). 

'^Calculated using B{F) = {N - Z) = 1 



the ^^C(g.s.)^-^'^N(g.s.)) might be due to A-isobar ad- 
mixtures in the reaction mechanism. 

In Rcfs. [17, 41], it was shown that the cross section 
measured to the 15.1 MeV state in the earlier experi- 
ments [13-15] was too high, possibly due to contamina- 
tion from the ^^C(p,n)^^N(g.s.) transition. Furthermore, 
it was shown in Ref. [41] that the cross sections mea- 
sured for ^•^C(p,n)^^N(15.1 MeV) reaction and its analog 
^^C(n,p)^'^B(g.s.) were consistent, providing further in- 
dication that the use of the former to calibrate the GT 
unit cross section is appropriate. 

Following Ref. [16], further arguments were provided 
in Ref. [17] that the proportionahty between GT strength 
and cross section for the ground-state transition was bro- 
ken, and that it should not be used to calibrate the GT 
unit cross section to extract B(GT) values for the excita- 
tion of other states in ^'^N. The authors of Ref. [17] quote 
a B(GT) for the transition to the excited state at 3.51 
MeV of 1.06±0.05, based on the empirical unit cross sec- 
tions obtained with even-even mass targets in the (p,n) 
reaction [15]. However, they also determined a GT unit 
cross section using the excitation to the 15.1 MeV state. 
Its value was 22% lower than that obtained using the em- 
pirical relationships. The reduced value for the GT unit 
cross section would increase the B(GT) for the transition 
to the excited state at 3.51 MeV to 1.29 ± 0.06, which 
is within error bars consistent to the value extracted via 
the (^IIe,t) reaction and close to the shell-model value. 

We should note that in both the {p,n) work discussed 
in Refs. [16, 17] and the {^He,t) analysis described here a 
large breaking of the proportionality between B(GT) and 
cross section for the ground state transition is reported. 



However, whereas in Refs. [16, 17] these are modeled 
through the modified operator in the shell-model calcu- 
lations mentioned above, we find that the breaking in 
the (''He,t) reaction at 420 MeV is a direct result of the 
tensor-T component in the effective nucleon-nucleon in- 
teraction. Since the consequences of either description 
are similar (causing interference between AL = and 
AL — 2 amplitudes) they could well represent the same 
underlying mechanism. 

We conclude that there is a significant spread in the 
B(GT) values reported for the transition to the 3.51 
MeV state in ^^N, but that the most recent results from 
^^C{p,n) and the new results from ^'^C('^He,i) favor an 
empirical value much closer to the shell-model calcula- 
tions than the earlier (p,n) results. Nevertheless, for the 
purpose of studying the effects on the pre-explosion evo- 
lutionary track of SNe la below, we used all available 
experimental values, as well as the shell-model results. 



V. ELECTRON-CAPTURE RATES IN THE SNe 
la ENVIRONMENT. 

As mentioned above, the electron capture on ^^N, pop- 
ulating states in ^"^C can be investigated using the analog 
transitions from ^'^C to ^'^N. We performed such a study. 
The logft for the ground state transition is known from 
/3-decay and thus fixed. Therefore, the only parameter 
that was varied in the study was the B(GT) for the tran- 
sition from the ground state of ^^N to the 3/2- state at 
3.68 MeV in -'^'^C. As mentioned above, isospin symmetry 
was assumed, so that the B(GT) values deduced from ex- 
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perimcnt using the transition from the ^"'C ground state 
to the 3/2" state at 3.51 MeV in could be applied 
directly. Electron captures onto ^''C are blocked by the 
high threshold for this reaction (13.4 MeV). 

The electron-capture rate calculations are based on 
the method described in Refs. [42 45], which were im- 
plemented into a new code [46]. The calculations were 
performed in a grid spanning logmipYe) values from 8.7 
to 9.6 in steps of 0.1 and dimensionless temperatures 
(Tg = T/109(K)) from 0.01 to 100. Electron chemical 
potentials {ne) were computed from a tabulation [47]. In 
Fig. 3a, the results of these calculations are shown for 
the case in which the B(GT) for the transition to the ex- 
cited state has been set to 1.5, which is the shell-model 
value, corrected for quenching as used in Ref. [5]. Near 
the temperature of relevance in the simmering stage of 
SNe la (approximately Tg ~ 0.4), the electron-capture 
rate is nearly independent of the temperature since kT 
is much smaller than the chemical potential He and thus 
only depends on the value of pYe- 

The relevant logjo(P^e) in the simmering stage ranges 
from about 8.8 to 9.5 and for the two extremes in this 
range, the electron-capture rates on '^'^N are plotted as a 
function of the B(GT) for the transition to the excited 
state at 3.68 MeV in Figs. 3b and c. The temperature 
Tg was fixed to 0.4 in both cases. If the B(GT)=0 for 
the transition to the 3.68 MeV state, capture to the ^^C 
ground state is the only possible channel. Prom Figs. 3b 
and c one can thus read off the branching ratio for cap- 
ture to the excited state as a function of its B(GT). At 
logiQ{pYe) = 8.8, the electron-capture rate is then 6.49 
s~^. If the B(GT) for the transition to the first excited 
state is taken from the older {p,n) data, the rate increases 
to 7.49±0.06 s~^. Using the newer {p,n) data, it increases 
to 7.78±0.06 s~^. The rates calculated based on the 
(^He,t) data and the shell-model calculation are close: 
8.15±0.09 s~^ and 8.31 s"^. Including the transition to 
the excited state increases the electron-capture rate, com- 
pared to that to the ground-state transition only, by 15%, 
20%, 26%, and 28% if the B(GT) values are taken from 
the older {p,n) data, the latest {p,n) data, the {^ile,t) 
data, and the shell-model calculation, respectively. 

At logio{pYe) = 9.5, the electron-capture rate has in- 
creased to 47.5 s~^ if the B(GT) of the transition to the 
excited state is set to zero. The contributions from the 
transition to the excited state increase that rate by an 
additional 30%, 39%, 51%, and 56% if the B(GT) val- 
ues arc taken from the older (p,n) data, the newer {p,n) 
data, the {^B.e,t) data, and the shell-model calculation, 
respectively. 

In a white dwarf star on the threshold of igniting and 
becoming a type la supernova, the timescale for a sound 
wave to traverse the star is thyd ~ {G{p))~^/'^ ~ Is, 
where G is the gravitational constant and {p) is the av- 
erage mass density of the white dwarf. The electron cap- 
ture rate is sufficiently fast, Aec < tj^yj, that during the 
pre-explosion "simmering phase," in which the heating 
from the fusion of ^^C-|-^^C occurs on a timescale > thyd, 




[ log,o(pYJ=8.7 

.2 .1 

10 10 1 10 
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0.25 0.5 0.75 1 1.25 1.5 1.75 2 
B(GT)(^^N(g.s.)-^(^^C(3.68 MeV)) 



FIG. 3: a) Calculated electron-capture rates on ^"'N in the 
stellar environment as a function of temperature and pYe 
(ranging from logio(pye)=8.7 to 9.6 in steps of 0.1). The 
B(GT) for the transition to the 3.68 MeV excited state in 
'^^C was fixed to the shell-model value of 1.5. b) Electron- 
capture rate at logio(pye)=8.8 and T9=0.4 as a function of 
the assumed B(GT) for the transition to the 3.68 MeV excited 
state in ^^C. The B(GT) values extracted from the (p,n) data 
("1" refers to Refs. [13-15] and "2" to Ref. [17]), (^He,t) data 
and the shell-model calculation are indicated and associated 
with electron-capture rates. The dashed lines correspond to 
error margins in the B(GT) values deduced experimentally, 
c) idem, but at logio(pye)=9.5. 



we should not expect the increased electron-capture rate 

to affect the reaction flows that occur during this sim- 
mering [5, 6]. As mentioned in the introduction, cap- 
tures into the excited state do affect, however, the heat 
deposited into the star from this reaction. Since the "'^^N 
is produced from ^^C(p, 7)^'^N with the protons coming 
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logio(pYe)=8.7 




0.25 0.5 0.75 1 1.25 1.5 1.75 2 
B(GT)(^^N(g.s.)-^(^^C(3.68 MeV)) 



FIG. 4: a) Calculated neutrino energy-loss rates due to the 
electron-capture on ^^N in the stellar environment as a func- 
tion of temperature and pYe (ranging from logio(pVe)=8.7 to 
9.6 in steps of 0.1). The B(GT) for the transition to the 3.68 
MeV excited state in ^^C was fixed to the shell-model value 
of 1.5. b) neutrino energy-loss rate at logio(pye)=8.8 and 
T9=0.4 as a function of the assumed B(GT) for the transi- 
tion to the 3.68 MeV excited state in ^'^C. The B(GT) values 
extracted from the (p,n) data ( "1" refers to Refs. [13-15] and 
"2" to Ref. [17]), (^He,t) data and the shell-model calcula- 
tion are indicated and associated with energy-loss rates. The 
dashed lines correspond to error margins in the B(GT) values 
deduced experimentally, c) idem, but at logio(pVe)=9.5. 



The heat evolved from, this reaction per second per ^^N 
nucleus is 

S = Aec {^ie + Q- E,) , (3) 

where the Q-value for the /'?+-decay is 2.22 MeV, and 
£^1/ is the mean energy of the neutrino. We compute 
El, = floss /Aec from the integration over the electron 
phase space, with Acc shown in Fig. 3 and i^ioss shown in 
Fig. 4. We find that, at logiQ{pYe) = 8.8, each reaction 
i3N(e-,i/e)^^C deposits into the white dwarf 0.53, 0.63, 
0.77, and 0.81 MeV for B(GT) taken from the older (p, n) 
data, the newer {p, n) data, the ("^He, t) data, and the 
shell-model calculation, respectively. At a higher den- 
sity, logio(pi^e) = 9.5, the heat deposition is 1.70, 1.89, 
2.08, and 2.12 MeV for these four cases. If the B(GT) 
were set to zero, the heat deposition from this reaction 
would be 0.07 and 0.99 MeV for log^o^pYe) = 8.8 and 9.5, 
respectively. 

We estimate the contribution of this reaction to the 
heating of the pre-explosive white dwarf by integrating 
a reaction network at logio(pY'e) = 8.8 and 9.5, with the 
temperature at a fiducial value of 0.4 GK. The integra- 
tion time was chosen to be the heating timescale CpT/S, 
where Cp is the specific heat and S is the heating rate. 
The dominant heat sources are the reactions ^^C + ^^C, 
i2C(n,7)i3c, i2c(p,7)13n and ^^G{a,ny^O, which re- 
lease about 2.7 MeV for each carbon consumed [5, 6]. At 
logio(P^e) = 8.8, the reaction ^'^N(e~, i^e)^^C contributes 
0.4%, 2.6%, 3.1%, 3.8%, and 4.0% for B(GT) = 0, from 
the older {p,n) data, the newer {p,n) data, from the 
(^He, t) data, and from the shell-model calculation, re- 
spectively. At logjQ(pFe) = 9.5, the higher electron Fermi 
energy ensures that this reaction contributes somewhat 
more to the heating, although the contribution is offset 
by a decrease in the production of ^^N. We find the con- 
tribution to the heating to be 4.3%, 7.0%, 7.8%, 8.5%, 
and 8.7% for B(GT) = 0, from the older {p, n) data, from 
the newer {p, n) data, from the (^He, t) data, and from 
the shell-model calculation, respectively. Although the 
differences between the results from the (^He, t) and the 
(p, n) data significantly affect the electron-capture rate 
and the net heat deposition in the star due to this tran- 
sition, we find that the overall effect on the pre-explosion 
evolution of the white dwarf is small. The capture rate 
is sufficiently fast in all cases that the rate of neutroniza- 
tion of the white dwarf is controlled by the production 
of "'^'^N, and at relevant densities the heating from this 
transition is a small component of the total. 



from the branch ^^C(^^C,p)^'^Na, the number of elec- 
tron captures onto ^"'N is proportional to the number 
of ^^C consumed in raising the white dwarf temperature 
to ignition. A greater heat release from ^^N(e~, z/e)^^C 
decreases the amount of carbon that must be consumed 
during this phase. 



VI. CONCLUSIONS 

In summary, we extracted the GT strength for the 
transition from ^"^C to the ^'^N 3/2~ excited state at 3.51 
MeV using the (^He,f) reaction at 420 MeV and existing 
data for the same reaction at 450 MeV. We find that the 
B(GT) is nearly consistent with shell-model calculations 
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using the CKII interaction and is higher than earUer ex- 
perimental values extracted using the ^'^G{p,n) reaction. 
However, the discrepancies between ('^He.f) and (j>,n) re- 
sults are lifted if, in the analysis of latest {p,n) experi- 
ments, the GT unit cross section is calibrated using the 
transition to the 3/2^ state at 15.1 MeV, as was done 
for the (^He,i) reaction. 

The Gamow- Teller component in the 
-'^^C('^He,t)^''N(g.s.) reaction was estimated to be 
strongly increased because of interference between 
AL = and AL = 2 amplitudes. As a result, the 
estimated Fermi component for this transition decreases. 
After this correction, both Fermi and GT unit cross 
sections extracted from the data, arc; consistent with 
the mass-dependent trends of Fermi and GT unit cross 
sections, respectively, measured in even-even nuclei. 

Under the assumption of isospin symmetry, the B(GT) 
for this transition is the same for the transition from 
the i3N(g.s.) to the 3/2" state at 3.68 MeV in ^^C. 
Combined with the known weak transition strength from 
the /3-decay of ^^N to the ground state of ^^C, we in- 
vestigated the impact of the differences between the ex- 
tracted B(GT) values from the (p,r?,) data, the ('^He,i) 
data and the shell-model calculation on the evolution of 



Type la supcrnovae during the prc-explosion simmering 
stage. Wc find that, although electron-capture and neu- 
trino energy-loss rates for this particular channel are sig- 
nificantly lower when the {p,n) values are adopted instead 
of the (^He,t) or shell-model values, the overall effect on 
the evolution is small, since the heating due to this re- 
action is small compared to that from other reactions. 
Nevertheless, the details of this reaction are now quanti- 
fied for use in further studies of SNe la. 
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